Mean-field phase diagram of cold lattice bosons in disordered potentials 
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We address the phase diagram of the disordered Bose-Hubbard model that has been reaUzed in several 
recent experiments in terms of optically trapped ultracold bosons. We show that a a systematic description 
of all of the expected quantum phases can be obtained at both zero and finite temperature via a site- 
dependent decoupling mean-field approach. Also, we relate the boundaries of the Mott insulating phase 
to an off-diagonal non-interacting Anderson model whose spectral features provide a new avenue for 
determining the debated nature of the phase surrounding the Mott lobes. Our approach is simple yet 
effective, and scalable to systems with experimentally relevant sizes and features. 
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PACS numbers: 03.75.Lm, 05.30.Jp, 64.60.Cn 

In 1998 Jaksch et al. demonstrated that the Bose- 
Hubbard (BH) model 10] could be accurately realized with 
a degenerate Bose gas in an optical lattice. Testament to 
the rapid pace of experimental development in this field, 
the defining superfiuid-insulator quantum phase transition 
of this model was first observed by Greiner et al. 13] in 
2002. The control over atoms afforded by optical lattices 
has served to initiate a broad range of investigations with 
bosonic and fermionic atoms |4, 5]. 

While optical lattices naturally provide a defect-free pe- 
riodic potential that allows precise control over the rela- 
tive strengths of interactions and tunneling, there is much 
interest relating to the interplay between interactions and 
disorder in the BH model. Indeed, experiments are cur- 
rently developing in this direction and have demonstrated 
several approaches for engineering disorder in optical lat- 
tices, such as: using laser speckle fields 0; an additional 
incommensurate lattice |7, 8, 9]; or using a distinguishable 
atom species to act as a randomly distributed set of impu- 
rities ll[i3]. 

As it was first demonstrated in Ref. |1], degenerate 
bosonic atoms in an optical lattice are described by the 
Bose-Hubbard (BH) Hamiltonian, 
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originally introduced as a toy model of superfluid ''He in 
porous media |2]. The operators n,„ = al^Qm, and 
al^ respectively count, destroy and create bosons at lattice 
site m, and obey canonical commutation rules [a™, a]^,] = 
^■m m' • The chemical potential /i is related to the total boson 
population = n^. The Hamiltonian parameters, 
namely the boson-boson (repulsive) interaction [/ > 0, 
and the hopping amplitude across neighbouring sites J are 



directly related to the atomic scattering length, and lattice 
depth |1]. The adjacency matrix Amm' appearing in the 
hopping term allows a simple algebraic description of the 
lattice structure, being finite if sites m and m' are nearest 
neighbours and zero otherwise. The local potential u,„ at 
site m is related to the features of the effective potential 
jfH]. Here this quantity will be random to realize disorder. 

As discussed in the seminal paper by Fisher et al. |2], 
the presence of disorder enriches the phase diagram of 
the BH model — that in the pure case (Vm = 0) con- 
sists of an extended superfluid region and a series of Mott- 
insulator lobes — with a further Bose-glass phase. Similar 
to the Mott insulator, the Bose glass phase is characterized 
by the absence of superfluidity, however has a finite com- 
pressibility (or gapless spectrum) like the superfluid phase. 
A representative list of techniques used to investigate the 
disordered BH pha se diagram includes field-theoretic ap- 
proaches Il2l lllill2l[l3|[ . decoupl ing (or Gutzwiller) mean- 
field approximatio ns jls j [141 llSl \l6i IitIi quantum Monte 
Carlo simulations JTMliaiand others llfl IH IH IH 24]. 
Nevertheless, several aspects of the problem are still sub- 
ject of active debate, such as a precise characterization of 
the different phases |24K the issue of the direct transition 
from MI to SF phase U Ellll, and the phase 

diagram at finite temperature iHtIi . 

In this paper we show that a site -dependent mean-field 
approach Ilia, captures all of the essential features of 
the phase diagram of the disordered BH model, at both zero 
and finite temperature. As summarized in the legend of fig- 
ure|4] the different phases, namely the Mott insulator (MI), 
the Bose glass (BG), the superfluid (SF) and — at finite 
temperature — the normal fluid (NF) are characterized by 
the value of three quantities, i.e. the superfluid fraction fs, 
the compressibility k and the condensate fraction f^. For 
simplicity we present results for a translationally invari- 
ant ID lattice with random on-site potentials Vm uniformly 
distributed in [—A, A]. However we emphasize that the 
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mean-field approach lends itself to more general situations, 
such as higher dimensional systems, different realizations 
of disorder and realistic trapping potentials ll2^ . Before 
describing our results we redefine the parameters using U 
as the energy scale, so that henceforth J, /x, f™, A are to 
be intended as J/U, fi/U, v^/U, A/U. 

The superfluid fraction is determined by the stiffness of 
the system under phase variations. 

Eg - Eq 

where (•) denotes thermal average and Eg is the energy 
of the system with twisted boundary conditions. The lat- 
ter are obtained by introducing the so called Peierls phase 
factors in the kinetic term of Hamiltonian Q. In ID this 
amounts to setting A„„„' = e^""" 5.„^„^,+i + e*^5,„,„'_i 
101 ■ The compressibility is defined as k = d^N = 
(3{{N'^) - {Nf), where (5 = U/keT is the inverse tem- 
perature. Finally, the condensate fraction, < /c < 1 is 
defined as the largest eigenvalue of the one body density 
')/A^ Hill- 



matrix (a.|„ a™,, J— I m. 

The decoupling mean-field approach IIM Hal results 
from the approximation aj^a^' ~ al^Om' + a'^cim' — 
ot*mOim', with a„i = {am)- This turns Hamiltonian (0 into 

^ = I]m=i "^m + h, where Hm = - Jhmal, + 
^^ara) and 
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Since the mean-field Hamiltonian Tl is the sum of on- 
site terms Tim, the original problem is reduced to a set 
of problems involving quantities relevant to neighbouring 
sites El El 



tr(a„ 



-I3H„ 



Indeed, according to Eq. @, the local Hamiltonian Tim de- 
pends on the a„j''s at sites m' adjacent to m, thus through 
the {a„i} spatial correlations in the system are approxi- 
mately included in the decoupling approach. One useful 
measure of these spatial correlations is the one-body den- 
sity matrix, = [{{rirn) - \am\^)dmm' + a*„^am']/N, 
as defined above. Note that the set of a^'s characterizing 
the state of the system can be seen as a stable fixed point 
of the map defined by Eq. l|5}. An easily found fixed point 
corresponds to the choice = for all sites. In Ref. 
the stability of such fixed point is studied for T > 0. In 
the T = case Eq. (JSjl turns into a,„ = (^|a„i|^') = 

{tpm\am\tpm), whcrc l^f) = 11™=! iV'm) and \tprn) are the 
ground-states of the entire system and of Tim, respectively. 
Note that in this limit the fixed point = corresponds 
to the number-squeezed ground state typical of the MI 
phase. Indeed, since Tim = Hm one easily gets \4'm) = 
\vm), where n,„|i/,„) = Vm\vm) and i/,„ = \^JL-Vm^ G N. 
Making use of first order perturbation theory it is possible 



to show that this MI phase is stable for J < J^ = \qM\ ^> 
where qM is the maximal eigenvalue of the matrix of en- 
tries Qmm' = gm{lj)^m,n', with (/i) = fll^ " ^m) 

and g{x) = (x + — x){x — \x\)\ ll2al . Note 

that for the a™ = fixed point fs = k = 0, while 
fc oc — > 0, as expected for a MI state (see leg- 
end. Fig. 0. As soon as am 7^ the local ground state 
is not a number state any more, \ipm) = Yl'^=o ^miyW) , 
and the system enters a compressible phase. Also, it can 
be shown that oc |amp, where the bar denotes average 
over the lattice sites. Hence the compressible phase found 
for J > Jc has a finite condensate fraction, i.e. long range 
correlations. Generally the superfluid fraction fs has to be 
evaluated numerically [E^. In the following we show that 
on ID systems the boundary of the SF region is simply re- 
lated to the vanishing of at some site of the lattice. Be- 
fore discussing our results, we rapidly review the pure BH 
model. Since = 0, Q = g{^)A, and the known equa- 
tion for the boundary of the Mott lobes, Jc = [2g{fj,)]^^ 
1 2], is easily recovered. As we mention above, k > as 
soon as a„ > 0. Furthermore, it is not hard to show that 
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V(^ni)' where the first equality is true in 



the thermodynamic limit M — s- oo. 

Let us now consider Vm uniformly distributed in 
[—A, A], focusing first of all on the boundaries of the MI 
phases. We begin by noting that the matrix Q whose max- 
imal eigenvalue gives the critical value J^, can be related 
to the Hamiltonian for an off-diagonal Anderson model 
whose random hopping amplitudes have an unusual dis- 
tribution [29] . This can be seen observing that Q has 
the same spectrum as the symmetric matrix of elements 

^mm' '^mm'^mm' \/ grn {f^)gm' (^t^)-^m m'; describ- 
ing noninteracting particles hopping across the sites of the 
lattice described by A with random amplitudes given by 
r. The spectrum of R can be analyzed for very large 
(M ~ 10^) ID systems using e.g. transfer matrix methods 
I30ll . It is easy to see that the evaluation of Jc = IqmI^ 
makes sense only for G J^, = + A, i' + 1 — A), 
where is a non negative integer and A < 1/2. Indeed if 
/i G = [i^ — A, + A] there is the possibility (certainty 
for M oo) that /i — = G N for one or more of 
the Vm's. This means that gmifJ-) diverges, and Jc = 0. 
Hence, as expected, the Mott lobes are found only within 
the ^ intervals J^, |2], where the entries of Q are always 
finite. We also note that if gmifj) is replaced by its average 
^(//), which discards the spatial correlations inherent in the 
Anderson model described by R our ap proach reproduces 
the boundaries given in Refs. H flTL 12111 . 

In what follows we discuss some numerical results ob- 
tained on ID lattices of size M = 100 where periodic 
boundary conditions are assumed. We consider a fixed 
realization for the "profile" of the disordered potential, 
Um G [—1,1] and obtain the actual potential at a given 
value of the strength A as Vm = An,„ . By averaging over 
disorder and considering large lattices we have verified that 
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FIG. 1 : Phase diagrams for a disordered lattice comprising M — 
100 sites for two values of the strength of the random potential, 
A = 0.05 (left) and A = 0.25 (right). 
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FIG. 2: Spectral density cr(e) of Q vs. ^ for a lattice compris- 
ing M = 10® sites. The leftmost (rightmost) panel corresponds 
A — 0.05 (A — 0.25). The spectra are normalized so that 
corresponds to e = 1. The small central panels show the corre- 
sponding phase diagrams, as indicated by the arrows. The gray 
areas are the SF regions, the solid black lines are the MI lobes. 



finite size effects are negligible in these results. The val- 
ues of /s have been obtained evaluating Eg with a self- 
consistent minimization procedure. We note that the same 
results can be derived from the self-consistent solution at 
= after a perturbative expansion in 9 l23l . Fig. ^shows 
the T = phase diagram of the system in the region of the 
first Mott lobe for A = 0.05 and A = 0.25. Note that 
the extended uniformly-colored region in the lower part of 
both panels refers to vanishing /g, according to the color- 
bar. The hatched regions correspond to the Mott lobes as 
evaluated computing the largest eigenvalue of the matrix 
Q. As expected, both k and obtained from the numerical 
minimization of the mean-field energy vanish inside and 
are finite outside the hatched regions. These results clearly 
show that the disordered potential induces the appearance 
of a phase that is absent in the pure model, characterized by 
finite compressibility (and condensate fraction) but vanish- 
ing superfluidity. This phase is hence naturally identified 
with the Bose-glass (BG). In more detail, increasing the 
strength of the disorder A causes the BG to extend at the 
expense of the SF and MI phases. In Ref. |2] Fisher et al. 
suggest that the MI to SF transition always occurs through 
a BG phase, a conjecture that has been actively debated, 
e.g. see llTa.[l3lfl^ ^1. The results presented here sug- 
gest that, in the mean-field picture, the MI and SF phases 
appear to be still connected near the tip of the Mott lobe for 
weak disorder (left panel of Fig. 0. Furthermore, our ap- 
proach provides a new avenue for understanding the pres- 
ence of a BG phase separating the MI and the SF via the 
spectral features of the Anderson model associated to the 
matrix Q. This can be seen in Fig. |21 showing the spec- 
tral density a of Q (or, equivalently, of R). Note indeed 
that for large disorder (right panel) a is always vanishing 
in the proximity of the band edge, whereas for small dis- 
order (left panel) one can recognize a clearly different be- 
haviour. In the vicinity of the tip of the lobe, where the SF 
and MI phases seem to be directly connected, a has an ev- 
ident peak, similar to what happens on a homogeneous ID 
lattice. We further observe that for ID systems /g vanishes 
as soon as the mean-field parameters vanish at two adjacent 
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FIG. 3: Square modulus (left) and phase differences (right) of 
Qfm' along the dashed line in the right panel of Fig^(A — 0.25, 
M = 0.35, e = 10"^). 



sites, e.g. ajn-i = = 0. Indeed in this situation one 
can verify that Eg = Eq since = |a^^|e*''", where 
the parenthetic superscript refers to the value of the Peierls 
phase and ip^n = '&m — "^m+i = Omodijn — m,M). 
This can be derived from the self-consistency constraint 



'a^Li + e'^a^^;^ observing that a^fj inherits 



the phase factor from 7^' = |7m^ |e"'" due to the specific 
form of Tim- Fig- |3] clearly illustrates this phenomenon 
displaying what happens on the lattice while crossing the 
dashed line in the right panel of Fig. Q for = 10^"^. As 
soon as |a^^ | (left panel) vanishes at two or more adjacent 
sites, /s = and ipm (right panel) stops fluctuating and 
equals 6 wherever it is defined. Of course on higher di- 
mensional systems the above argument does not apply, and 
one expects the onset of superfluidity to be related to the 
percolation of the a„i's through the lattice 1 16]. 

The last results we present here are the same phase di- 
agram as in Fig. ^ but at a finite temperature T = 
O.OlU/kB- As it can be seen in Fig. |4] the values of 
fs, K and fc allow the characterization of four different 
phases. Strictly speaking, at finite temperatures the MI is 
replaced by a normal fluid (NF) which is always compress- 
ible. However at small temperatures some regions of the 
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FIG. 4: The same phase diagrams as in Fig^ at a finite tem- 
perature T = 0.01[//fcB. Filled and empty circles in the legend 
denote finite and zero values for the relevant quantity. 



NF phase feature a compressibility so small that they can 
be considered (quasi) MI 1 14, 25, 28]. In particular our re- 
sults show that at finite T the BG cannot be characterized 
simply as a compressible non-superfiuid phase. Indeed this 
is true also of the NF phase, where /s = because a,„ = 
everywhere due to thermal fluctuations, even in the absence 
of disorder. The condensate fraction distinguishes between 
the NF and BG, i.e. the region = is comprised of 
BG phase with /c 7^ — where the superfluidity is de- 
stroyed (mainly) by the presence of disorder — and the NF 
phase with /c = - where the superfluidity is destroyed 
(mainly) by thermal fluctuations. As T is increased the 
NF dominates over the BG phase, while the (quasi) incom- 
pressible MI lobes shrink, until eventually only SF and NF 
phases remain. We mention that the latter phase has glassy 
features according to the treatment described in Ref. UtIi . 

In summary, in this work we show that the site- 
dependent mean-field decoupling approach captures the 
phases of the disordered Bose-Hubbard model at both zero 
and finite temperature, and that in general three indicators 
should be taken into account to characterize the phase di- 
agram. We observe that the boundaries of the MI lobes 
can be related to a non-interacting off-diagonal Anderson 
model. In particular, we present the phase diagrams of 
ideal one-dimensional systems for different values of the 
strength of the disorder, and suggest that the specific na- 
ture of the transition from the MI to the SF phase might 
be related to the spectral features of the relevant Ander- 
son model, which we investigate for lattice sizes up to 
M = 10®. Also, we observe that on ID systems the 
superfluidity is destroyed as soon as the local mean-field 
parameters vanish at two adjacent sites, and sketch an ar- 
gument explaining this. We conclude observing that the 
site-dependent decoupling mean-field approach appears to 
be a very promising tool for investigating more general 
situations, including different realizations of the disorder, 
higher dimensions and realistic features, such as the har- 
monic confinement that characterizes actual experiments 

SEil]. 
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